We present a search for low-mass brown dwarfs in the 
Pleiades open cluster. The identification of Pleiades members 
fainter and cooler than those currently known allows us to 
constrain evolutionary models for L dwarfs and to extend the 
study of the cluster mass function to lower masses. We con- 
ducted a 1.8 deg 2 near-infrared /-band survey at the 3.5 m 
Calar Alto Telescope, with completeness J cp i ~ 19.0. The de- 
tected sources were correlated with those of previously avail- 
able optical /-band images (I cp i ~ 22). Using a J versus I - J 
colour-magnitude diagram, we identified 18 faint red L-type 
candidates, with magnitudes 17.4 < J < 19.7 and colours 
I - J > 3.2. If Pleiades members, their masses would span 
~0.040-0.020 Mq. We performed follow-up HK s -b<md imag- 
■ ing to further confirm their cluster membership by photometry 
and proper motion. Out of 1 1 1 J candidates with proper motion 
^vq • measurements, we find six cluster members, two non-members 
\ and three whose membership is uncertain and depends on the 
^ intrinsic velocity dispersion of Pleiades brown dwarfs. This 
dispersion (>4 mas yr -1 ) is at least four times that of cluster 
stars with masses ^1 Mq. Five of the seven other IJ candi- 
dates are discarded because their J — K s colours are bluer than 
those of confirmed members. Our least massive proper motion 
| members are BRB 28 and 29 (~25 Mj up ). The J versus I - J 
J> sequence of the L-type candidates at J > 18 is not as red as 
| theoretical models predict; it rather follows the field L-dwarf 
sequence translated to the cluster distance. This sequence over- 
lapping, also observed in the J versus J - H and J - K di- 
agrams, suggests that Pleiades and field L dwarfs may have 
similar spectral energy distributions and luminosities, and thus 
possibly similar radii. Also, we find a — 0.5 + 0.2 for a power- 
law approximation dN/dM oc M~ a of the survey mass spec- 
trum in the mass range 0.5-0.026 Mq. This value is similar 
to that of much younger clusters, indicating no significant dif- 
ferential evaporation of low-mass Pleiades members relative to 
more massive ones. 
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Abstract. 



1. Introduction 

As a result of progressive cooling, brown dwarfs are expected 
to evolve from late M spectral type at very early stages (age 
$ 10 Myr) to L- type (r eff in the range -2400 to -1400 K; 
iDahn et alj E()02)/ and at sufficiently old ages (a few Gyr), to 
T-type and beyond (r e ff $ 1 400 K). The mass-spectral type 
and mass-luminosity relationships at different ages remain to 
be derived from observations of cool dwarfs, which also per- 
mit the c alibration of the theory of substellar evolution (e.g. 
IChabrier' & Baraffel EoOOl) . In comparison to the spectra of M 
dwarfs, L dwarfs are characterized in the optical by the weak- 
ening of the metal-oxide TiO and VO bands due to the con- 
densation of the metals Ti and V in dust grains that may sig- 
nificantly affect the atmospheric structur e and the emergen t 
spectrum (e.g. the DUSTY models from lChabrier et all200ol) . 
Bands of metal-hydrides CrH, FeH and CaH and bands of wa- 
ter become more intense especially in the near-infrared (IR), 
whereas the neutral alk ali metals Na, K, Rb, Cs and Li becom e 
stronger in the optical (Chabrier & Baraffe 2000; Basri 2000). 
The emergence in the near-IR of CH4 bands characterizes the 
beginning of the T dwarf sequence. Field L dwarfs have I - J 
and J - K s colours redder than ~3 . 3 and ~1 .3 mag, respect ively, 
incre asing from early to later spectral subclasses ( Martin et"al] 
Il999h . For the L4-L8 spectral subclasses the near-IR J — H 
and J - Kg colours app e ar to saturate at ~1.2 and ~2.0 mag 
jKirkpatrick et alJ l2000t iBurgasser et ai1 Eo02). respectively, 
whereas for T dwarfs these colours become blu er, with val- 
ues decreasing to below zero (K napp et all 2004 ). From L0 to 
T8 spectral types the I - J colour is increasing to ~5.8 mag 
JPahn et alJ2002h . Associated with the other colours it can give 
a rough indication of the spectral type. 

The Pleiades open cluster (~120 Myr, ~130 pc) served 
for decades as a reference stellar laboratory where models 
were contrasted with observations. Since the discovery of 



Pleiades brown dwarfs llRebolo et alJl 19951 [l 996 ), subsequent 
studies have provided numerous fainter substellar candidates 
with spectral types down to late M llZapatero Osorio et al. 
1 1 997 M iBonvier et al l ll 998l hereafter B98: iFesfinl 11998? 
IZapatero Osorio et alJ 1 19991). from which a fraction is con- 
firmed by lithium detection ( Stauffer et al. 1998b; M artin et alJ 
l2000h or by proper motion jMoraux et alJ l200l[ hereafter 
M0 1). Very large area proper m otion surveys jHamblv et alJ 
1999; Deacon & Hamblv 2004) have obtained on the other 
hand a significant census of the stellar-substellar boundary pop- 
ulation of the cluster. The coolest object for which a spec- 
tral type h as been obtained is the L0 brown dwarf candidate 
Roque 25 IMartrn et alJll998bl) . with an estimated theoretical 
mass of ~0.035 M0. A dozen fainte r objects than Roque 25 
have been identified in deep surveys ( Festin 1998a; Beiar 2000; 
iDobbie et alJ2002UNagashima et all2003UMoraux et alJ2 003 , 
hereafter M03; Schwartz & Becklin 2005), and some were con- 
firmed by proper motion (Bo uv et alJ2 006). 

The mass spectrum dN/dM = f(M) of the Pleiades 
cluster, where dN stands for the number of objects in the 
mass range dM, is best f itted by a lognormal function jM03t 
Dea con & Hamblvll2004l) . and its approximation by a power 
law M~ a in the low mass range 0.6-0.03 Mp provid es a spec- 
tral index a ~ 0.6-0.8 (Dobbi e et all 120021 lM03). Similar 
values are found in other much younger open clusters, as 
for example cr Orionis (~2- 4 Myr), where a — 0.8 + 0.4 
for 0.2 < M(M Q ) < 0.013 feeiar et al.ll200ll) . and IC 348 
(~3 Myr) wher e a = 0.7 + 0.2 for 0.5 < M(M ) < 0.035 
(Tei et al. 2002). It is important to determine the behaviour of 
the a index in the lower mass domain because it can help to 
discriminate among brown dwarf formation mechanisms. For 
instance, for the formati on of grav itationall y unstab le cores by 
turbulent fragmentation JPadoan & Nordlundll2004l) . the mass 
spectrum is predicted to be log-normal, with a peak that ap- 
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pears at smaller masses when the sonic Mach number and the 
mean density of the cores are greater. 

In this paper we present a 1.8 deg 2 deep /-band survey 
of the Pleiades open cluster aimed to find lower-mass brown 
dwarfs (<0.035 M@), b y com parison with a new analysis of 
/-band data obtained by B98. Our goal is to identify L dwarf 
cluster candidates and study the location of the L dwarf se- 
quence in the colour-magnitude diagram. Preliminary results 
have been reported in Biha in et"afl d2005l) . Because the cluster 
has a well known age, the properties of Pleiades L-type brown 
dwarfs could provide useful constraints to models predicting 
the evolution of substellar objects. Additionally, we can esti- 
mate the cluster mass spectrum of the surveyed area at lower 
masses. 

2. Observations and data reduction 

2.1. J -band survey 

We obtained /-band images on 1998 October 27-28 with the 
Omega-Prime (£2') instrument at the 3.5 m Telescope of the 
Centra Astronomico Hispano Aleman (CAHA; Calar Alto, 
Spain). The fields were specifically chosen to cover the opti- 
cal Rl-band survey from B98 obtained with the 3.6 m Canada- 
France-Hawaii Telescope (CFHT)/UH8K survey. We observed 
the fields CFHT 2, 5 (north half), 6, 7, 8, 10, 11 and 13 (see 
their spatial distribution in Fig. [3- They correspond to an 
area of the cluster of ~1.8 square degrees. The Q' camera has 
a 1024x1024 HgCdTe Hawaii detector with a pixel scale of 
0.40" pixeL 1 and a field of view of 6.8' x 6.8'. To cover each of 
the 29'x29' UH8K fields (hereafter CFHT fields, as designed in 
lB98h a connecting observation macro was used. The exposure 
times and repetitions in each pointing were 16 s x 2, except for 
CFHT7 and CFHT8 (16 s x 3), and CFHT6 (18 s x 2). Because 
the offset between successive macro pointings corresponded to 
one third of an Q.' field of view, each area was repeated in to- 
tal 6-9 times (96-144 s), depending of the CFHT field. The 
individual images were sky-subtracted, flat-fielded using dome 
flats, combined and aligned in horizontal strips of ~ 30' x 6.8', 
using standard scripts and routines within the IRAF 1 environ- 
ment. These were later analyzed for aperture and PSF photom- 
etry using DAOPHOT. The PSF photometry was calibrated with 
measurements from the Two Micron All Sky Survey (2MASS) 
Point Source Catalog dCutri et alJl2003h and which had errors 
of less than 0. 1 mag. On average, about five 2MASS calibrators 
were used per horizontal strip. The /-band completeness (lim- 
iting) magnitude was estimated at the maximum (fainter half 
maximum) of the histogram of object counts per magnitude. 
We found values of ~19.0 and ~19.5 mag, respe ctivel y. 

The optical survey in the RI bands from B98J was ob- 
tained with the UH8K instrument at the CFHT (Mauna Kea, 
Hawai'i) on 1996 December 9-13. The UH8K camera com- 
prises eight Loral 2048x4096 CCD detectors with a pixel scale 
of 0.21" pixel . These data were used initially for a study of 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. Location of the surveyed Pleiades fields (squares) rela- 
tive to the brightest Pleiades members with visual magnitude 
V < 5 (large open circles). The small open circles represent the 
brown dwarf candidates from B98 with lithium or proper mo- 
tion consistent with cluster membership, and the filled circles 
represent all the other candidates identified in this study. 



brown dwarfs and low-mass stars by B98. In the presen t pa- 
per we use an independent data analysis by Beiar (2000). The 
photometric calibration of the CFHT fields was performed by 
comparing our photometry to the Cousins photometry of their 
candidates and applying an average offset; for the fields CFHT6 
and CFHT8 however, different offsets were applied because 
they were observed through thin cirrus. The completeness and 
limiting magnitudes were ~23.5 and -25.0 mag in the R band, 
and ~22.5 and -23.5 mag in the / band, respectively. 

To i dent ify fainter an d redder objects than those found 
by EH and lBeiarl |2000), we used as reference the detected 
/-band sources and searched for their optical counterparts in 
the /-band images using a FORTRAN correlation program, 
CORREL, kindly provided by M. R. Zapatero Osorio. The cor- 
related objects were plotted in a / versus / - / diagram (see 
Fig. |3 and contrast ed wit h the 120 Myr DUSTY isochrone 
from Chabrieret al. (2000), shifted to the cluster distance of 
133.8+3.0 pc JPercival et alJl2005l) . As most of our faint ob- 
jects were bluer than the theoretical isochrone in the / ver- 
sus / - / diagram, we decided to establish our selection of 
candidates based on an empirical criterion. We considered a 
bluer envelope to cluster candidates (dashed line in Fig.|2ji, de- 
fined by the blu est lith ium or prop er motion Pleiades objects 
reported JMartm et alJl200(j iMOll) : CFHT-P1-1, Calar 1 and 
CFHT-P1-25 (horizontal thick marks in Fig. 13 . The envelope 
was extended to faint magnitudes / > 18.5 at a fixed colour 
of / - / = 3.3, to avoid most of the M8-9 dwarf contami- 
nants which have bluer colours. As potential cluster members 
we have selected 34 objects (named "BRB") with / > 14.4 
and colours 1.7 < / - / < 3.8 (Table EJ. We checked that 
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Table 1. RIJ photometry and coordinates for the low-mass stars and brown dwarf candidates. 



Name 


Names reported'' 


R-I± 0- fi _; 


1 ±(Tj 


J ±(Tj 


/ - J 


KA (JzUUU) Dec. {Sl\)\)\)) 

{ I- „\ to t ft \ 

(Ami) ( ) 


BRB 1 


CFHT-P1-1, BPL 242, MBSC 91, (*) 


1.76 ± 0.13 


16.10 ± 0.07 


14.41 ± 0.03" 


1.69 


03 51 51.6 +23 34 50.2 


BRB 2 


CFHT-P1-2, BPL 267, DH 765, (*) 


2.00 ± 0.10 


16.57 ± 0.07 


14.66 ± 0.06 


1.91 


03 52 44.3 +23 54 15.2 


BRB 3 


CFHT-P1-6, (*) 


2.38 ± 0.10 


17.13 ± 0.07 


14.71 ± 0.05 


2.41 


03 52 07.9 +23 59 14.6 


BRB 4 


(*) 


2.38 ± 0.10 


17.03 ± 0.07 


14.73 ± 0.04 


2.30 


03 44 23.2 +25 38 44.7 


BRB 5 


CFHT-P1-3, HHJ 22, BPL 272, MBSC 99, (*) 


1.94 ± 0.10 


16.66 ± 0.07 


14.81 ± 0.06 


1 .85 


03 52 51.8 +23 33 48.9 


BRB 6 


CFHT-P1-4, BPL 280, MBSC 101, (*) 


1.92 ± 0.10 


16.85 ± 0.07 


14.96 ± 0.06 


1.89 


03 53 09.6 +23 33 48.3 


BRB 7 


CFHT-P1-5, DH 590, (*) 


2.11 ± 0.10 


16.99 ± 0.07 


14.98 ± 0.06 


2.01 


03 48 44.7 +24 37 22.7 


BRB 8 


CFHT-P1-7, BPL 253, MBSC 108, (*) 


1.84 ± 0.13 


17.48 ± 0.07 


15.16 ± 0.07 


2.32 


03 52 05.8 +24 17 31.7 


BRB 9 


CFHT-PM2, BPL 294, CFHT-PLIZ-6, (*) 


2.55 ± 0.10 


18.00 ± 0.07 


15.20 ± 0.03 


2.80 


03 53 55.1 +23 23 37.4 


BRB 10 


CFHT-P1-9, BPL 202, MHOBD 6, (*) 


2.17 ± 0.10 


17.78 ± 0.07 


15.46 ± 0.07 


2.32 


03 49 15.1 +24 36 22.4 


BRB 11 


CFHT-P1-13, Teide 2, BPL 254, CFHT-PLIZ-3, (*) 


2.11 ± 0.13 


18.14 ±0.07 


15.50 ± 0.07 


2.64 


03 52 06.7 +24 16 01.4 


BRB 12 


CFHT-PH 1, Roque 16, BPL 152, (*) 


2.25 ± 0.10 


17.92 ± 0.07 


15.64 ± 0.07 


2.28 


03 47 39.0 +24 36 22.1 


BRB 13 


CFHT-P1-15, (*) 


2.44 ± 0.10 


18.65 ± 0.07 


15.99 ± 0.07 


2.66 


03 55 12.5 +23 17 38.0 


BRB 14 


CFHT-P1-21, Calar 3, BPL 235, CFHT-PLIZ-12, (*) 


2.49 ±0.11 


18.93 ± 0.07 


16.13 ± 0.05 


2.80 


03 51 25.6 +23 45 20.6 


BRB 15 


CFHT-P1-25, BPL 303, CFHT-PLIZ-20, (*) 


2.70 ± 0.13 


19.68 ± 0.07 


16.71 ± 0.04 


2.98 


03 54 05.3 +23 34 00.2 


BRB 16 


PIZ 1 


— 


20.10 ± 0.12' 


16.84 ± 0.07 


3.26 


03 48 31.4 +24 34 37.7 


BRB 17" 


(**) 


- 


20.92 ± 0.07 


17.42 ± 0.06 


3.50 


03 54 08.31 +23 54 33.4 


BRB 18" 


CFHT-PLIZ-28, (**) 


2.37 ± 0.17 


21.20 ± 0.07 


17.61 ± 0.07 


3.59 


03 54 14.08 +23 17 52.2 


BRB 19" 


— 


2.25 ± 0.14 


20.95 ± 0.09 


17.79 ± 0.05 


3.16 


03 54 51.49 +23 45 12.2 


BRB 20" 


CFHT-PLIZ-35, (*), (**) 


2.55 ± 0.23 


21.47 ± 0.06 


18.06 ± 0.07 


3.40 


03 52 39.16 +24 46 29.7 


BRB 21" 


(**) 


2.34 ± 0.21 


21.68 ± 0.08 


18.14 ± 0.05 


3.54 


03 54 10.27 +23 41 40.3 


BRB 22" 


CFHT-PLIZ 2141 


2.14 ± 0.20 


21.97 ± 0.08 


18.31 ± 0.05 


3.65 


03 44 31.27 +25 35 15.1 


BRB 23" 


(**) 


2.13 ± 0.27 


22.03 ± 0.10 


18.55 ± 0.10 


3.49 


03 50 39.53 +25 02 54.5 


BRB 24 


(*) 


2.95 ± 0.51 


22.01 ± 0.09 


18.71 ± 0.08 


3.30 


03 48 19.65 +24 52 09.5 


BRB 25 






22.24 ± 0.09 


18.74 ± 0.07 


3.50 


03 52 59.73 +23 51 56.0 


BRB 26 




1 .yj ± U.OU 


00 AO j- A AQ 


lis, /o ± u.iu 


i 1 1 

.1 . 5 1 


AO AQ 17 AO /10 1 
UJ +o JO. 1 / +ZO Uy +0. 1 


BRB 27" 


CFHT-PLIZ 1262, (**) 


2.30 ± 0.32 


22.64 ±0.11 


18.89 ±0.09 


3.74 


03 44 27.24+25 44 41.9 


BRB 28" 


(*), (**) 


2.64 ± 0.45 


22.36 ± 0.08 


19.02 ±0.10 


3.33 


03 52 54.92 +24 37 18.6 


BRB 29" 






22.84 ±0.10 


19.05 ± 0.07 


3.79 


03 54 01.43 +23 49 58.1 


BRB 30 




1.98 ±0.23 


22.48 ±0.10 


19.17 ±0.07 


3.31 


03 47 56.96 +24 28 58.4 


BRB 31 




2.50 ± 0.25 


23.05 ±0.10 


19.29 ±0.10 


3.76 


03 51 25.93 +24 15 32.2 


BRB 32 






22.72 ±0.12 


19.36 ± 0.09 


3.36 


03 49 00.86 +24 41 38.5 


BRB 33" 






23.31 ±0.20 


19.63 ± 0.09 


3.68 


03 51 20.15 +23 45 18.4 


BRB 34 






23.01 ±0.16 


19.64 ± 0.09 


3.37 


03 48 46.55 +24 45 03.2 



" With follow-up observatio n for p roper motion measurement (see Tablel2*l. 

b (*) IBeiai200a). CFHT-P1 1593). BPL IPinfield et aU200CD. MBSC and CFH T-PLIZ |M0l. (**) [Bihain et all2005l). DH ID eacon & Ham bhl2004tl . CFHT-PLIZ IBouy et alj 
1200a HHJ 22 pamblv et al!1993fl . MHOBD 6 lstaufrerSll998ill . Teide 2 IMarti'n et alll998at) . Rogue 16 Izapatero Osorio et a7ll997rir Calar 3 IMartin et 31119961) . PIZ 1 

Cossburn et al. 1997). In Fig.fTlwe overplotted the 14 CFHT-P1 objects (ope n circles) and the 20 other objects (filled circles). 
c Coordinates of CFH T-P1, CFHT-PLIZ, PI Z 1 and the other objects are from B98 , M03 , C ossburn et aT] Q99JI and the present study, respectively. 
d /-band estimate from Martin et al. ' 2000 ) because our estimate is bad pixel-contaminated. 

e /-band estimate corrected for a ~ 0.25 mag reddening due to a close star and ray light traces ; the offset was obtained by comparing to objects of similar counts in clean sky areas of the 
image. Converting the Kitt Peak filter measurement lap = 19.64 from Cossburn et al. ' 1997) to Cousins with the transform equation given by Jameson et al. ' 2002). we find / = 19.9. 



these objects had colours R — I k, 1.8. Some fainter candidates 
had unreliable or too faint R-bmd magnitudes and no value is 
given in Tabled We discarded as contaminants the previously 
identified non-p roper m otion candidates CFHT-P1-8 (|M01) and 
CFHT-PLIZ-25 dM03h . 

Of the 34 BRB candida tes, 22 we re identified already by 
previous surveys (e.g. Beiar 2000, B98; see footnote of Tabled 
for all the references). In the magnitude range of the brightest 
objects we identify an additional candidate, BRB 4, which is 
also in the 2MASS catalog. This object is not mentioned by 
B98 and is not in the list of proper motion and Rl -b and pho- 
tometric candidates from D eacon & Hamblvl J2004I) . The lat- 
ter survey overlaps the CFHT survey and its /?-band limiting 
magnitude is greater than the magnitude of BRB 4, R = 19.4. 
Therefore this object should have been detected. Possibly it was 
blended or it is not a proper motion cluster member. 

The magnitude errors of the BRB candidates are 0.06 mag 
for J ~ 14.5 - 17.0 and 0.08 mag for J ~ 17.0 - 20.0. 
Comparin g our 7-band magnitudes with those available in the 
literature (M artin et aDl2000t Ipiiifield et alJl2003l) . i.e. those in 
the range J ~ 14.5 - 17.0, we find: 

1- < -/brb - -/Martin > = -0.04 + 0.09, for the objects in com- 
mon CFHT-P1-2, 3, 4, 5, 6, 7, 9, 11, 12, 15, 21 and 25. 



CFHT-P1-6 and CFHT-P1-7 present the greatest differences, 
-0. 18 and -0.26 mag, respec tively. Both objects lack of Hot 
emission (Mar tin et alJl200oT) . but CFHT-P1-7 has a proper 
motion consistent with that of the cluster. In one of our 
individual /-band images, the radial profile of CFHT-P1- 
7 peaks at relatively higher counts and is not well cen- 
tred, probably due to a cosmic ray; this might explain the 
greater brightness (it is also supported by the difference 
with 2MASS 7 RR r - 7 2 mass = -0.14). For CFHT-P1-6, as 
discussed by M01, if this is an equal mass binary as sug- 
gested by its position well above the cluster sequence in the 
colour-magnitude diagram (B98; see also Fig. the bina- 
rity might affect the determination of its short term proper 
motion and therefore explain its deviance from the cluster 
peculiar motion. The binarity could also explain some of 
the magnitude difference observed. 
2. < /brb - /pinfieid > = 0.04 ± 0.09, for CFHT-P1-2, 3, 4, 
7, 9, 11, 12, 13, 21 and 25. CFHT-P1-4 and CFHT-P1-11 
present the greatest differences, 0. 19 mag and 0. 15 mag, re- 
spectively. These are also found for 2MASS: 7brb - Jimass 
= 0.11 and 0.19, respectively. But CFHT-P1-4 and CFHT- 
Pl-11 have no deviant proper motions or anomalous posi- 
tions in the colour-magnitude diagram. 
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Fig. 2. / versus I- J colour-magnitude diagram for the correlated survey objects. The circular symbols correspond to 77 candidate 
or confirmed cluster members (Tabled, with redder I - J colours than the dashed line boundary (see text for details). At / 
> 17.4 we plot: proper motion members (large filled circles), possible proper motion members (circled filled circles) and non- 
prope r motion m embers (large empty circles). The solid line and the dash-dotted line represent the -120 Myr DUSTY isochrone 
( Chab rier et alJ [2000) and the ~125 Myr NextGen isochrone (Baraffe et al. 1998), respectively, shifted to the distance of the 
cluster. Masses in solar units are indicated for the DUSTY isochrone. The triangles correspond to field dwarfs shifted to the 
distance of the cluster. Finally, the upper and lower dotted lines at the bottom of the diagram indicate the completeness and 
limiting magnitudes of the survey, respectively. 
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Our /-band magnitudes rely on an approximate offset to 
the magnitudes from B98 and have an error of ~0.1 mag. 
The /-band magnitude from B98 is obtained after transforma- 
tion onto the standard Cousins system, with observation of red 
Landolt standards. Due to the lack of very red standard stars, 
the transformation is linearly extrapolated for the red and faint 
objects and could produce systematic errors. C ompar ing our 
magnitudes to those in the Cousins system from M03, for the 
candidates in common and 1-2 mag fainter than CFHT-P1-25 
(one of the faintest objects from B98), we observe a difference 
I ere - ^Moraux ~ 0.1, which is quite small in comparison with 
the I - J > 3.3 threshold that we applied to select the faint 
candidates. 

Finally, from the 34 BRB Pleiades candidates, the 18 
faintest have / ^ 17 and / — / S 3.3, the magnitudes and 
colours expected for L-type dwarfs. 

2.2. Follow-up observations 

2.2.1. /T s -band imaging 

Follow-up /if s -band imaging with CAIN-II at the 1.55 m 
Telescopio Carlos Sanchez (TCS; Teide Observatory, Tenerife) 
was performed for all 1 8 //-band photometric substellar candi- 
dates (see the observation log in Table[2ll. The CAIN-II camera 
has a 256x256 HgCdTe Nicmos3 detector with a pixel scale of 
1 .00" pixeL 1 in wide field configuration, providing a 4.2' X4.2' 
field of view. The observation macro consisted in 6 s exposures 
X 10 repetitions x 10 dither positions. Total exposure times 
spanned a range of 10-90 min, depending on the expected 
brightness of the candidates and the seeing. Data were sky- 
subtracted, flat-fielded (using sky flats instead of dome flats), 
aligned and combined in a similar way as for the /-band data 
(Sect. 12. It . We also performed the photometry and estimated 
the completeness and limiting magnitudes as described for the 
/-band. On average, about six 2MASS calibrators were used 
per field. Average completeness and limiting magnitudes dur- 
ing the different nights are indicated in Table|2 

2.2.2. HK s -band imaging for astrometry 

For the proper motion determination of the candidates, we ob- 
tained subarcsecond H- and /T s -band images with the Long- 
slit Intermediate Resolution Infrared Spectrograph (LIRIS) at 
the 4.2 m William Hershel Telescope (WHT; Roque de los 
Muchachos Observatory, La Palma) and //-band images with 
Q2000 at the 3.5 m Telescope of CAHA (see Tabled). 

We d ecided also to obs erve the Pleiades brown dwarf 
Teide 1 jRebolo et all 1 1995b with both LIRIS and 02000 
to improve its proper motion measurement and to check if 
there are any systematic errors in the proper motion measure- 
ments using different instruments. Our first epoch image was a 
12'x 12' reduced image, obtained with the TEK3 CCD detec- 
tor (0.7" pixeL 1 ) at the prime focus of the 2.5 m Isaac Newton 
Telescope (INT; Roque de los Muchachos Observatory, La 
Palma) on 1995 December 19 (Zapatero Osorio 2005, private 
communication). 



The LIRIS camera has a 1024x1024 HgCdTe Hawaii de- 
tector with a pixel scale of 0.25" pixeL 1 and a 4.2' x 4.2' field 
of view, whereas the Q2000 camera has a 2048x2048 HgCdTe 
Hawaii2 detector with a pixel scale of 0.45" pixeL 1 and a field 
of view of 15.4' x 15.4'. The observation macros consisted of: 
5 s exposures x 18 repetitions x 5 dithers (LIRIS night 2005 
January 23), 5 s exposures x 12 repetitions x 5 dithers (LIRIS; 
January 24), 3 s x 20 coadds x 7 dithers (Q2000; February 
1), 20 s x 6 repetitions x 5 dithers (LIRIS; March 25), 2 s 
x 20 coadds x 15 dithers (02000; October 21, 22 and 23), 
2 s x 15 coadds x 20 dithers (Q2000; October 24) and 3 s 
x 15 coadds x 20 dithers (Q2000; October 25). The total ex- 
posure times were chosen to achieve enough S/N (^20 in the 
peak) for precise measurement of the position of the candi- 
dates. They ranged between 5 and 40 minutes, depending on 
the expected brightness of the candidates and the seeing. 

The raw images were first bad-pixel corrected using a bad- 
pixel mask and the IRAF routine PR0T0.FIXPIX. The mask 
was obtained from the flat (sky-flat for LIRIS and dome-flat 
for £22000) using N0A0 . IMRED . CCDRED . CCDMASK. The im- 
ages were then sky-subtracted - sky images were obtained by 
combining 5-20 consecutive images of similar sky counts - 
and divided by the flat. 

The resulting images were distortion corrected, aligned us- 
ing 10-20 reference stars and combined. For LIRIS these last 
three steps were performed with the LDEDITHER task within 
the package IRAF . LIRIS . L IMAGE developed by J. A. Acosta- 
Pulido. For Q2000, the images were distortion corrected 
by projecting them on the celestial coordinate grid with 
MSCIMAGE. The grid was obtained by astrometry with the 
USNO-A2 catalog and the script myasrtrom.cl from E. Puddu. 
This script invokes CCFIND, that correlates the catalog sources 
with those in the image, and CCMAP, which computes the 
plate solution. Usually, 50-100 stars were accepted by a 
Legendre order-4 transformation fit and a rejection threshold of 
1 .5 sigma, providing a projection transformation with standard 
deviation ~0.12" in both right ascension and declination. The 
distortion corrected images of £22000 were then aligned and 
combined using the INTALIGN task within the LIRIS package. 

The Z//f s -band photometry was obtained similarly as for 
the /-band data. On average, six and eleven 2MASS calibra- 
tors were used per LIRIS and Q2000 field, respectively. The 
average completeness and limiting magnitudes during the dif- 
ferent nights are indicated in Tabled 

2.3. Proper motion analysis 

The proper motion measurements were obtained as follows. 
First we selected non-saturated and well defined (S/N > 10) 
single objects among all those that appeared within ~3' of 
a brown dwarf candidate. We measured their centres in the 
first- and the second-epoch images with the IRAF task CENTER 
(DIGIPH0T.DA0PH0T package), with a precision of 0.01- 
0.1 pixels in both image dimensions. For BRB 22 and 33 we 
used the //-band £22000 images as second-epoch images, be- 
cause they provided greater S/N for the candidates than the 
£22000 Zf s -band images. Then we computed the transformation 
from second-epoch image to first-epoch image with GE0MAP. 
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Table 2. HK s -band follow-up observations. 



Telescope/instrument 


Filter(s) 


Date 


ant/n't 

UDject^s ) 


< c>, < / >" 


Tf- Q IC* A TXT TT 


is 

a s 


ZUU4 


uec / 


D D D 1 "7 10 Ol TO OA 

riKri 1 /, lo, zl, zj, zu 


1 1 fx 
~1 /.O, 


1 C 
1 o.Z 


TPQ /r A TM TT 
1 L-o/L-AlfN — 11 


V 
* S 


ZUU4 




dKd ZU, Z5 


1 B 1 
~15.1, 


~~ 1 O.O 


WHl/LlKlvS 


J_J 

ti 


ZUUj 


Jan 2 j 


ODD n 11 

riKri 1 /, Zj 


~ iv.j, 


on n 
-ZU.U 


\\ ' I IT /I TDTC 

Wnl/LlKlo 


ti 


ZUUj 


T„„ 1/1 

Jan 24 


D D D OA OQ 

riKri ZU, zo 


1 n n 


^ iy.j 


J.j m Calar/iizUUU 


ti 


ZUUj 


reu i 


D D D 1 "7 1 Q Ol T Q « ^ a. 1 

dKd 1/, lo, zl, leiae 1 


1 c ^ 

~ lo.J, 


1 Q T 


Tr^C If - A TM TT 

1 L.O/L.A1JN — 11 




ZUUj 


Mar 6 


D D D 

riKri zj 


1 f ^ 


1 7 n 
v l / .U 


\\ 71 IT /I TDTC 

WHI/LIKLS 




ZUUj 


Mar zj 


D D D 0"7 Tal^a 1 

dKd z / , leiae 1 


10 

~ lo.o, 


■^iy.j 


Tf^C Z^" 1 A TXT TT 

1 L-o/LAlfN — 11 


is 

» s 


ZUUj 


n„» i t 

yjci i / 


D D d in 
riKri iy 


1 1 ^ 

~1 /.J, 


1 c ^ 
~~ lo.J 


Tf^C tf~* A TXT TT 

1 L.O/L.A1JN — 11 


is 

£ s 


ZUUj 


UCt 15 


D D D TI 

riKri LL 


1 c n 
lo.U, 


1 o n 
-~ iv.u 


j.j m L,alar/i,2zUUU 


ti 


mn^ 
ZUUj 


Uct 21 


D D d in 
riKri IV 


1 n n 




Tf^C /f" 1 A TXT TT 


A s 


zuuj 


\JCl Zl 


D D D 1A 1A 

riKri Z4, j4 


~ 1 / .0, 


lo.o 


j.j m L,aiai/i,2zuuu 


XT, A s 


2005 


Oct 22 


Dnn 

rSKrS ZZ 


-19.0, 


-19.5 


TCS/CAIN - II 


*S 


2005 


Oct 22 


BRB 29 


-18.0, 


-19.0 


3.5 m Calar/n2000 


H 


2005 


Oct 23 


BRB 29 


-18.5, 


-19.0 


TCS/CAIN - II 


K s 


2005 


Oct 23 


BRB 31, 33 


-18.8, 


-19.5 


3.5 m Calar/f22000 


K s 


2005 


Oct 24 


BRB 33 


-19.0, 


-19.5 


3.5 m Calar/fi2000 


H 


2005 


Oct 25 


BRB 33 


-19.0, 


-19.5 


TCS/CAIN - II 


K s 


2005 


Oct 28 


BRB 30, 32 


-18.0, 


-19.0 



" Average completeness and limiting magnitudes. 



For the Q2000 image of Teide 1 , we performed a transforma- 
tion from first- to second-epoch (i.e. from TEK3 to Q2000), 
because the destination image had a wider field and a smaller 
pixel scale and therefore provided more accurate equatorial as- 
trometry with the USNO-A2 catalog. With GE0MAP we selected 
the objects whose positions minimized the transformation error 
and could thus be considered as reference objects. A Legendre 
function of order 3 was used when more than 12 objects were 
available, and of order 2 for fewer objects. The error, related 
to the scatter of these objects, ranged over 0.03-0.15 pixels. 
For the candidate BRB 23 the error in the y dimension was 
as high as 0.6 pixels, because of vertical smearing affecting 
the object positions in the /-band image (field CFHT 6). With 
GEOXYTRAN we used the transformation to predict the candi- 
date position (x\p,yi p ) in a first-epoch image, and compared 
it with the measured one, (xi,y\), obtaining the pixel shifts 
(A Xl ,A Vl ) = (x\ p - X\,y\ p — y\). The plate solution of a first- 
epoch (or destination) image was obtained with the USNO-A2 
catalog and the myasrtrom.cl script, in a similar way as ex- 
plained above. Then using WCSCTRAN we converted the shifts in 
pixels to shifts in equatorial coordinates. Due to the precision of 
the equatorial astrometry (<x - 0.12"), its error contribution to 
the proper motion was negligible. But for the LIRIS proper mo- 
tion measurement of Teide 1 , the best first epoch TEK3 fit had 
a relatively high standard deviation (ct R a = 0.6", o~ dec = 0.7") 
and the corresponding proper motion value differed by almost 
~1 mas yr~' with those derived from less precise fits; thus we 
added -1 mas yr~' to the proper motion error in this case. In 
general, we assumed that the proper motion error was : 

I ""—* V ») = { Ai ' At j ' (1) 

where cr vl vl x2 , yi are the errors of the (x,y) pixel positions in 
the first- (1) and second-epoch (2) images, o- lx ly the error of the 
transformation between these images, 

,cost5, ii S the transforma- 
tion from (x, y) to (RA, DEC) coordinates, and At the time in 
Julian years between the two epochs of observation. 



3. Proper motion Pleiades low-mass brown dwarfs 

3.1. Photometry and proper motion results 

The HK S magnitudes for the red and faint brown dwarf candi- 
dates that we have followed-up are given in Table [3] The H- 
band magnitude of BRB 17 is the error-weighted average from 
the values obtained with £22000 and LIRIS, and the ^ s -band 
magnitudes of BRB 22 and 33 are the error-weighted averages 
from the values obtained with O2000 and TCS. Our 18 follow- 
up objects have magnitudes 20.9 < I < 23.3, 17.4 < J < 19.6, 
and colours 3.2 < / - J < 3.8, 0.7 < J - H < 1.3 and 
0.9 < J - K s < 2.2. Most of them remain good L dwarf 
candidates, what supports our I - J selection criterion. But 
BRB 24, 26, 30, 3 1 and 34 present much bluer J - K s colours 
than cluster proper motion members with similar /-band mag- 
nitude (see Fig. |3- Moreover, if they were late L or T dwarfs, 
they would certainly present much redder I - J colours (>4- 
5 mag). We will assume that these are non-cluster members. 
Finally BRB 32, with its red J - K s and relatively blue I - J 
colours, remains a possible cluster member. 

The proper motion results obtained for 11 of the brown 
dwarf candidates and two lithium brown dwarfs (Teide 1 
and Calar 3) are listed in Table |4] and are represented by 
small cir cles in the vector point diagram of Figure El Because 
Calar 3 (Ma rtin et all 1 19961 iRebolo et all 1 19961) is at a few 
arcminutes from BRB 33, we could also measure its proper 
motion and found a very good agreement with the mea- 
surement of ImOIL We note that the double-check measure- 
ments of Teide 1 and BRB 17 are fully consistent within 
the statistic errors. In Fig. [3] are represented also: (/) the ref- 
erence objects used to estimate the proper motions (aster- 
isks), (if) the circle of radius three times the greatest compo- 
nent of the average proper motion error derived from Table |4] 
((< (r^coss >, < cr m >) = (3.0, 4.3) mas yr~'). This circle 
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Table 3. Follow-up HK S photometry of the low-mass brown dwarf candidates. 



Name 


J 


I - J ± <Tl-j 


H ± cr H 


J - H ± (Tj-h 


K s ± cr Ks 


J - K s ± <Tj_ Ks 


BRB 17° 


17.42 


3.50+0.09 


16.75+0.02 


0.67+0.06 


16.20+0.06 


1.22+0.08 


BRB 18" 


17.61 


3.59±0.10 


16.74+0.06 


0.87+0.09 


16.08+0.07 


1.53+0.09 


BRB 19" 


17.79 


3.16±0.10 


16.99+0.05 


0.80+0.07 


16.67+0.05 


1.12+0.07 


BRB 20" 


18.06 


3.40+0.09 


17.33+0.05 


0.73+0.08 


16.56+0.07 


1.50+0.10 


BRB 21" 


18.14 


3.54±0.09 


17.05+0.07 


1.09+0.09 


16.39+0.07 


1.75+0.09 


BRB 22" 


18.31 


3.65±0.09 


17.38+0.04 


0.94+0.06 


16.69+0.04 


1.62+0.06 


BRB 23" 


18.55 


3.49±0.14 


17.47+0.06 


1.08+0.11 


16.65+0.06 


1.90+0.11 


BRB 24 


18.71 


3.30±0.12 






17.49+0.07 


1.22+0.11 


BRB 25 


18.74 


3.50±0.12 






16.91+0.24 


1.83+0.25 


BRB 26 


18.76 


3.31+0.13 






17.81+0.10 


0.95+0.14 


BRB 27" 


18.89 


3.74±0.14 






17.16+0.09 


1.74+0.12 


BRB 28" 


19.02 


3.33+0.12 


17.89+0.05 


1.13+0.11 


17.00+0.08 


2.02+0.13 


BRB 29" 


19.05 


3.79+0.12 


17.76+0.05 


1.29+0.09 


16.88+0.08 


2.17+0.11 


BRB 30 


19.17 


3.31+0.12 






17.65+0.09 


1.52+0.11 


BRB 31 


19.29 


3.76+0.14 






18.20+0.15 


1.09+0.18 


BRB 32 


19.36 


3.36+0.15 






17.42+0.20 


1.94+0.22 


BRB 33" 


19.63 


3.68+0.22 


18.71+0.07 


0.92+0.11 


17.87+0.06 


1.76+0.11 


BRB 34 


19.64 


3.37+0.18 






18.8+0.5 


0.9+0.5 



" With proper motion measurement. 



Table 4. Proper motions of the low-mass brown dwarf candi- 
dates and the lithium brown dwarfs Teide 1 and Calar 3. 



Name 


Instrument 


p„cos<5 ± cT^osa 








(mas yr"" 1 ) 


(mas yr~') 


Teide 1 


LIRIS 


22.5+3.6 


-40.5+4.0 


Teide 1 


£22000 


19.2+2.5 


-45.5+3.4 


Calar 3 


£22000 


24.5+2.2 


-32.6+2.2 


BRB 17 


LIRIS 


20.9+1.5 


-22.8+2.4 


BRB 17 


£22000 


17.7+2.9 


-20.7+3.1 


BRB 18 


£22000 


28.4+2.7 


-37.4+2.7 


BRB 19 


£22000 


75.7+1.0 


-143.9+7.7 


BRB 20 


LIRIS 


21.8+3.0 


-48.1+2.0 


BRB 21 


£22000 


18.1+3.4 


-49.1+4.1 


BRB 22 


£22000 


18.8+3.5 


-51.9+2.9 


BRB 23 


LIRIS 


29.2+3.0 


-25.6+14.7 


BRB 27 


LIRIS 


35.3+3.7 


-48.8+2.0 


BRB 28 


LIRIS 


15.1+3.1 


-37.4+2.1 


BRB 29 


£22000 


26.4+3.3 


-48.8+2.9 


BRB 33 


£22000 


-4.1+3.0 


-5.2+3.7 



is centred on the cluster average proper motio n, Gu„cosd>, ut) 
= (19 .15 + 0.23, -45.72 + 0.18) mas yr 1 (Robi chon et all 
Il999h . and we will assume that the objects within are cluster 
members. 

The L dwarf candidates BRB 17, 23 and 27 lie outside the 
circle in Fig. [3] and their membership is considered uncertain. 
The proper motion dispersion of the seven objects within this 
circle is (5.2, 5.9) mas yr~' . Subtracted quadratically by the av- 
erage proper motion error, it yields (4.3, 4.1) mas yr -1 , as an 
estimate 2 of the intrinsic velocity dispersion of Pleiades brown 
dwarfs. This velocity dispersion is at leas t four times that of 
Pleiades stars with masses £1 Mq (LTonesJll970t lpi nfield et alJ 

2 At the cluster distance, it corresponds to a tangential velocity of 
-2.5 kms 1 . 



Il998t) and appears consistent with the linear relationship of 
equipartition of energy, between the velocity dispersion and 
the inverse square root of the mass for cluster members , as 
expected in an nearly relaxed cluster llPinfieldetalJ [l998 and 
Fig. 4 therein). 

We cannot rule out that some of the IJ candidates lo- 
cated close to the region surrounding the circle in Fig. [5] are 
in fact cluster members (although field contaminants are nei- 
ther excluded, see Sect. 14.11 and Table [3}. This is actually 
the case of the Pleiades lithium brown dwarf Calar 3, sug- 
gesting that the intrinsic velocity dispersion of brown dwarfs 
could be even greater than our previous estimate. We need 
further data (for instance from lithium observations) to study 
the membership of th e pec uliar IJ candidates. Note also the 
case of CFHT-P1-15 (|B98), which has a proper motion of 
(66.5 + 8.1, -54.1 + 8.1) mas yr 1 ( iMOll) . thus located far 
from the cluster point in the proper motion diagram, although it 
presents a radial velocity consistent with that of the cluster, Li 
absorption, H„ emission a nd a red I - K s colour JStauffer et alJ 
ll998bHMartm et all2000l) . 

However, BRB 19 and BRB 33 appear clearly to be non- 
proper motion members; BRB 19 has an especially high 
proper motion. Finally, our error-weighted average estimate for 
Teide 1 is (20.3 + 2.1, -43.4 + 2.6) mas yr -1 , very close to 
the cluster p roper motion. It also agrees with the measurement 
obtained by iRebolo et al.1 dl995h . the only previous estimate 
(16 + 9, -35 + 16) mas yr , but with large conservative errors. 
Thes e measur ements, together with the clear lithium detection 
dRebolo et alJIl 996h . thus fully confirm Teide 1 as a Pleiades 
brown dwarf. Also, the photometry obtained with £22000 and 
LIRIS, H = 15.54 + 0.06 and = 15.0 7 + 0.09, respectively, 
agrees with that of ljameson et al.ld2002l) . H = 15.65 + 0.09 and 
K= 15.08 + 0.05. 
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Fig. 3. Vector point diagram of proper motion for the candi- 
dates (small circles with BRB number indicated) and the ref- 
erence objects (asterisks). The results for the Pleiades lithium 
brown dwarfs Teide 1 and Calar 3 are also shown. Double- 
check objects are linked by a solid line. The plotted large cir- 
cle is centred on the Pleiades average proper motion (see text 
for details). The reference objects of BRB 23 (small squares) 
present a large dispersion in declination which causes the large 
error bar of this candidate. The arrow plotted for BRB 19 indi- 
cates a location outside the diagram. 



14 r 



15 



16 



17 - 



18 - 



19 



20 t_ 
0.0 




Fig. 4. J versus J - H colour-magnitude diagram for all 
the Pleiades very low-mass star and brown dwarf candi- 
dates with lithium or proper motion consistent with mem- 
bership (filled circles). Their photometry was compiled from 
the present study. iMartm et alJ J2000t). th e 2 MASS catalog 
Zana tero Osorio et alJJl997bl) . lBeiarl (l2000) and ljameson et alJ 
The circular symbols (J > 17.4), lines and triangles are 
defined as in Figure |2 The typical error bars for the probable 
and possible cluster members by proper motion are represented 
on the left. 



3.2. Pleiades L dwarf sequence 

In the J versus I- J colour-magnitude diagram (Fig.|2j, the six 
cluster proper motion members (large filled circles) define the 
Pleiades L dwarf sequence. Note that the three possible mem- 
bers with relatively high proper motion (circled filled circles) 
share the same sequence (as well as in the near-IR colour- 
magnitude diagrams, Figs.|4]and|5}. 

In Fig. 13 we also present theoretical predictions for the / 
and J bandpasses. The solid line corresponds to the DUSTY 
isochrone and th e dash-dotted fi ne cor responds to the NextGen 
isochrone from Bar affe et alJ fl998) for an age 125 Myr. 
Masses in solar units are shown for the first isochrone. The 
NextGen model provides a good fit for J < 16.0-16.5 and the 
DUSTY model provides a good fit for 7^18. Both models fail 
to reproduce the photometric sequence at fainter magnitudes; 
the former does not account for the onset of dust formation in 
cooler atmospheres, whereas the latter overestimates the I — J 
colour. 



Also in Fig. |2 we show the sequence of field dwarfs (tri- 
angles) with parallactic distances, translated to the distance of 
the Pleiades cluster. Th ese nearby cool dwarfs are compiled 
from the literature by Caballero (2006). Their /- and JHK S - 
band photometry are from iDahn et alJ d2000l 2002) and the 
2MASS Point Source Catalo g, respectivel y. For the M dwarfs, 
the parallaxes are from Dahn et al. (2000, 2002) and the spec- 
tral types ar e fr om iGo limowsk i et al.1 (120041) . ILeggett et alJ 
(200ol l2002l) and lDahn et alJ d2002t). For t he L a n d T dwarfs 
the pa rall axes are from |Perrvman et al.l (1997), Dahnetal. 
d2002|) andlVrba et all J2004I) , and th e spectral type s are m ostly 
from IVrba et all d2004 . else from ILeggett et afl J2002h and 
Geb alle et al.l (1200 11) . In Fig. |2]the spectral types range over 
~M3-T2 from top to bottom. At magnitudes J < 17 the field 
dwarfs are bluer than the cluster brown dwarfs, but at J > 17 
the field sequence begins to intercept that of the cluster and to 
follow it slightly towards redder colours. 

In the J versus J - H and J - K diagrams (Figs.|4]and[5}, 
we show the Pleiades very low-mass stars and brown dwarf 
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Fig. 5. / versus J-K colour-magnitude diagram, with the same 
definition of symbols and error bars as in Figure |4] The small 
empty circles represent the IJ candidates without proper mo- 
tion measurements. 



candidates having lithium in their atmospheres or proper mo- 
tions consistent with membership in the cluster (filled circles). 
The circular symbols (/ > 17.4), lines and triangles are defined 
as in Figure |5] The DUSTY and NextGen isochrones provide 
relatively good colour predictions for J > 17 and J < 17, re- 
spectively, in both diagrams. The fact that the DUSTY model 
agrees with the observations in the / versus J - H and J-K 
diagrams but not in the J versus I- J diagram supports the pos- 
sible underestimation of the far red flux (/ band) relative to the 
near-IR flux (J-, H- a nd K bands). In Fig, g] the M9 Pleiades 
brown dwarf NPL 40 (Festin 1998a b) has a very blue J - H 
colour, as noted by Pinfield et al. (2003). In Fig. [5] the brown 
dwarf CFHT-P1-7, with its magnitude J = 15.16 probably en- 
hanced by a cosmic ray, is positioned at the left of the upper 
sequence, and the five very blue candidates at / > 18.5 (empty 
circles) are probable contaminants. 

In these near-IR colour-magnitude diagrams, where field 
dwarfs are also represented (~M3-L6), similar sequence over- 
laps appear at J > 17 as in the J versus I - J diagram. Thus 
the sequence of the L brown dwarf candidates of the Pleiades 
is nearly indistinguishable from that of the field dwarfs shifted 
to the cluster distance. Note that the average errors in the dis- 
tances of the sample field dwarfs are < <x d > = 0.4 pc (M5-M9) 
and 1 .0 pc (L0-L7), and the average errors in their translated J- 
band magnitudes, including the photometric errors and the dis- 



tance error of 3 pc for the Pleiades cluster, is of only 0.08 and 
0. 12 mag, respectively. The average error in their I - J colours 
is 0.03 mag (M5-M9) and 0.06 mag (L0-L7). These magni- 
tude and colour errors are small in comparison to the (field and 
cluster) trends observed. 

From the J versus I - J, J - H and J-K diagrams it ap- 
pears that, opposite to the Pleiades late M dwarfs, many of 
the Pleiades L dwarf candidates have the same colours and 
absolute magnitudes as their field counterparts. This suggests 
that they may have the same spectral energy distributions and 
luminosities. If the effective temperatures estimated from the 
spectral energy distribution do not depend much on the gravity 
(as indicated by spectral synthesis), then, they may have also 
the same effective temper atures. These obje cts would therefore 
have the same radii. From lAllen et al.l d2005h - the mean age of 
field L dwarfs is in the range 3-4 Gyr and the mean mass is in 
the range 0.06-0.07 Mq, which according to the DUSTY mod- 
els imply radii of ~0.09 Rq. These objects have almost reached 
the end of the contraction. We note however, that according 
to the same models (and using the luminosity estimates from 
Sect. the radii of the Pleiades L brown dwarf candidates 
are of ~0.125 Rq, a value inconsistent with the previous one. 
Either the predicted evolution of radii with time may have to be 
revisited by the models, so the Pleiades L dwarfs reach radii of 
0.09 Rq in about 120 Myr or, the estimated masses and ages of 
the field L dwarfs may require a revision toward higher masses 
or smaller ages. A determination of the radii from the lumi- 
nosities and the effective temperatures require precise spectro- 
scopic observations which are beyond the scope of the present 
study. 

4. Substellar luminosity function 

The /-band luminosity function, the number of objects versus 
bin of /-band magnitude, is represented in Fig.|6]for BRB 1- 
34, our Pleiades very low-mass star and brown dwarf //-band 
candidates (/ < 18: solid line, / > 18: dotted line). The error 
bars are Poissonian errors, i.e. the square root of the number of 
candidates per bin. The vertical dotted line delimits the com- 
pleteness of the survey, at / ~ 18.75. For this and fainter mag- 
nitudes, the present luminosity function is incomplete (dotted 
line). 

4.1. Contaminants 

The number of field M and L dwarfs that may contaminate our 
sample is obtained as follows. First we subdivide the region in 
the / versus /- / diagram occupied by the 34 objects of Tabled 
in small rectangular sections. These sections have heights of 
half /-band magnitude and widths delimited by the bluer enve- 
lope (dashed line in Fig.fJJ and the reddest object. 

Then, we consider the typical /-/colours of M- and L-type 
dwarfs as a function of the spectral type, based on the same 
sample of nearby cool dwarfs as mentioned in Section I3.2I 
Contaminants in our s urvey are es sentially M5-L4 dwarfs. 
Using the relation from lDahn et alJ J2002) between the spec- 
tral type of M7-L8 dwarfs and the maximum distance at which 
they can be observed for a given /-band magnitude (a relation 
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Fig. 6. Luminosity function for all the survey candidates (solid 
and dotted lines) and after contaminant correction (dashed 
line). The vertical dotted line represents the completeness limit 
of the survey. 



CFHT-PLIZ-25 and CFHT-PLIZ-25) and in Col. "Confirmed", 
the BRB numbers of the proper motion and lithium Pleiades 
members. In the next three columns, we give the estimated 
numbers of contaminants for the spectral type ranges M5-M7, 
M8-M9 and L0-L4. In column "total", the total number of con- 
taminants is provided for the entire spectral type range M5- 
L4. Nine contaminants are alr eady found (Col. "Con taminants 
found"): (i) CFHT-P1-8 dMOlh . CFHT-PLIZ-25 (M03). BRB 19 
and BRB 33, (ii) BRB 24, 26, 30 31 and 34. The first set of ob- 
jects is found by proper motion measurements and the second 
set by / - K s photometry (Col. "Measurement"). The numbers 
of contaminants found in each spectral type range are indicated 
in parenthesis. Subtracting these from the total number of con- 
taminants, we obtain the number of remaining contaminants 
that we may find among our survey candidates (Col. "iV r . c ."). 

After correcting the number of 77 candidates for contami- 
nants in each /-mag range, we find an agreement with the num- 
ber of confirmed Pleiades members, except for 15.0-16.0 and 
19.0-20.0 mag. In these ranges one more contaminant than ob- 
served is predicted. For the latter range, this can be explained 
by its incompleteness. In both cases, we assume the observed 
number of confirmed members as the effective value when de- 
riving the effective luminosity and mass functions, thus de- 
creasing by one the value of N r . c .. In Fig. |6] the effective lu- 
minosity function is represented by the dashed histogram. 



that we assume also valid for M5 and M6 dwarfs), we com- 
pute for any spectral type contaminating a given rectangular 
section the pair of distances associated to the /-band magni- 
tude boundaries of the section. Each pair of distances defines 
a volume towards the Pleiades, e.g. a truncated right pyramid 
with a square base, with top and base surfaces subtended by 
the survey angular area, 1.8 deg 2 . This volume is located at the 
average of the pair of distances (< d >). 

We use the densiti es of dw a rfs in the solar neighbourhood 
provided by Kirknatrick et al. (1994) for M5 and M6 spec- 
tral types, and bv lCruz et al.1 (bOOaT for M7-L4 at d < 20 pc. 
We take into account the completeness and sky coverage of 
these surveys. IChen et all J200ll) studied the star counts of 
the Sloan Digital Sky Survey for several hundred thousands 
stars at high latitudes north and south of the Galactic plane. 
They obtain an exponential scale height of the old thin disk of 
330+ 3 pc (late-type stars) and a Sun's distance from the plane 
of 27+ 4 pc. We assume that this scale height applies also for 
early L dwarfs. To compute the numerical density of a given 
spectral type dwarf at an average distance < d > from the Sun, 
we multiply the solar neighbourhood density by an exponential 
exp-(| < d > sin(fopi) -z©|/zhX where zh is the scale height, 
Zq the Sun's distance to the Galactic plane, bp\ = -23.52° the 
Pleiades cluster's galactic latitude. 

Finally we multiply the volumes by the corresponding den- 
sities and obtain the statistical numbers of contaminants. Each 
section in the / versus I - J diagram is thus contaminated by 
field dwarfs of different spectral types at different distances and 
with different densities. In Table|5] we summarize our estimate 
of contaminants. In the first column we indicate the /-mag 
range, in Col. "iV can ", the number of IJ candidates (including 



5. Substellar mass function 

Finally, we present the implications of our results in the sub- 
stellar mass function of the cluster surveyed area. The bolo- 
metric luminosity and effective temperature predictions of the 
DUSTY and NextGen mod els are mostly i ndependent of the 
atmospheric properties dChabrier et all200oT) . in contrast to the 
predicted magnitudes (e.g. for the M-L transition), which are 
not. We decided to infer a bolometric luminosity using a /-band 
bolometric correction and then to compare it with the model 
prediction for the distance and age of the Pleiades to obtain the 
mass. 

Ad opting the Pleiades cluster distance from lPercival et all 
d2005l) . we convert the /-band magnitudes to absolute mag- 
nitudes My, and add a bolometric correction BCj depending 
on / - /. This bolometric correction is obtained bv lCaba llero 
(2006J) from a fit to (M bo i - Mj, I - /) of the nearby cool 
dwarfs mentioned in Section 13.21 Binaries, peculiar dwarfs 
and dwarfs with line emission are not taken into account in 
the fit. For the nearby M dwarfs, the bolometric magnitudes 
are fromlGolimowski et all (bOOAlLeggett et alJ J20001 120021) 
and iDahn et alJ J2002h. For the near by L and T dwarfs, these 
are m ostly from IVrba et alJ d2004l) . else from iLeggett et alJ 
d2002l) and IGeballe et all J200ll) . The bolometric magnitudes 
Mboi = Mj + BCj that we obtain ed are then normalized to the 
solar value M bol o = 4.74 (ILiving ston 2002) and converted to 
luminosities. Using these luminosities, we interpolate linearly 
the theoretical data points (DUSTY model for BRB 8-34 and 
NextGen model for BRB 1-7) and obtain the masses. These 
are grouped into five bins: two low-mass stellar and three sub- 
stellar bins (see Table |6}. The lowest mass bin is incomplete 
(/ ^ 18.75). The differences between the masses derived from 
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Table 5. Contaminants among our low-mass star and substellar candidates. 



J range 


Wen" 


Confirmed 6 
(Li or pm) 


Estimated numbers of contaminants 
M5-M7 M8-M9 L0-L4 total 


Contaminants 
found 


Measurement 


N t .c. c 


14.5-15.0 


l d 


1, 2, 3, 5, 6, 7 


0.3 






0.3 









15.0-16.0 


1 


8, 9, 10, 11, 12, 13 


1.6(1) 






1.7 


CFHT-P1-18 


prop. mot. 


1 


16.0-17.0 


3 


14, 15 


0.8 


0.3 


0.1 


1.2 






1 


17.0-18.0 


4 


17, 18 




1.3(2) 


0.3 


1.6 


CFHT-PLIZ-25,BRB 19 


prop. mot. 





18.0-19.0 


8 


20, 21, 22, 23, 27 




1.7(2) 


1.3 


3.0 


BRB 24, 26 


J-K, 


1 


(19.0-20.0) 


7 


28, 29 




2.0(2) 


3.8(2) 


5.9 


BRB 30, 31, 33, 34 


J - K s , prop. mot. 


2 



" Including CFHT-P1-18 and CFHT-PLIZ-25. 
* BRB number. 

c We assume only and 1 contaminants that remain to be found in the J ranges 15.0-16.0 and 19.0-20.0, respectively, when deriving the 
effective luminosity and mass functions (see text, Sect. 14. 1> . 
d Including BRB 1 which has J = 14.41. 



Table 6. Mass bins for our mass function. 



Ay 


AM 


A'brb 


Wont ^eff 




(Mq) 






14.41-14.97 


0.147-0.104 


6 


6 


14.97-15.57 


0.104-0.075 


5 


5 


15.57-17.51 


0.075-0.039 


6 


1 5 


17.51-18.75 


0.039-0.026 


8 


3 5 


(18.75-19.65) 


0.026-0.018 


9 


6 3 



the bolometric luminosities of the NextGen- and DUSTY mod- 
els for the fainter objects BRB 8-22, which extend into the 
lower luminosity range of the NextGen model, are smaller than 
0.002 Mq and do not affect the distribution of objects in these 
bins. Each bin contains numbers of BRB candidates (A^brb) 
and contaminants (N CO nt), identifiable from Tables^and[5] The 
width of the mass bins within the completeness of the survey 
is such that each bin contains approximately the same effective 
number of objects, N e g = /Vbrb - N cont - 

In Fig.0 we show our results in the mass spectrum repre- 
sentation, log dNI dM, log M. Error bars refer to the Poissonian 
error. We show also the bins for the ~50 massive proper mo- 



, r 



in 



tion Pleiades objects from Deacon & Hamblv (2004) which are 
present in our fields 3 . For the ~35 fainter objects, correspond- 
ing essentially to M dwarfs (T e ff £ 3600 K), we obtain the 
masses with the method explain ed above, but using /- and J- 
band magnitudes from Deacon & Hamblv ( 2004) and 2MASS, 
respectively, because these objects are saturated in our images. 
For the brighter ones we proceed similarly but, this time, us- 
ing a synthetic Zf-band bolometric c orrection, function of the 
(J - K)cn colour and determined by Houda shelt et alJ d2000l) 
for FGK stars with 4000 K < r eff < 6500 K. We assume 
solar metallicity and log g — 4. The CIT colour is transformed 
into the 2MASS equivalent (see the Explanatory Supplement to 
the 2MASS All Sky Data Release 4 , Cutri et al. 2006) and the 
bolometric correction is applied on 2MASS 7/T s -band photom- 
etry. For the objects in the K-M transition, the masses obtained 




3 The objects in common with our sample, DH 765 and DH 590, 
are taken into account only once and in our mass spectrum bins. 

4 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.html 



-1.0 -0.5 
logM(M sol ) 

Fig. 7. The derived mass function for the Pleiades survey area. 
The filled circles are for the BRB obje cts, whereas the squares 
are for the objects from iDeacon & Hamblvl J2004I) . The right 
dotted line indicates the hydrogen-burning mass limit. The 
lower limit arrow for the lowest mass bin indicates that this 
bin is incomple te. The curve represents th e log normal func- 
tion obtained bv lDeacon & Hamblvl (|2004) scaled down to the 
total number of objects in the survey area and between 0.5 Mq 
and the completeness mass limit 0.026 Mq (56). The dashed 
line represents the power law fit to the data points between 0.5 
and 0.026 Mq, with a = 0.5 + 0.2. 



by any of these two methods are within the highest mass bin 
(M > 0.5 Mq), thus fixing the number of objects in this bin. 
A linear fit to all the data points between 0.5 and 0.026 Mq 
provides an a — 0.5 + 0.2 for the power law dN/dM oc M~ a . 
It is compatible with the mass bin beyond the completeness 
limit of the survey (lower limit arrow at the left of the vertical 
dotted line, Fig. [7}, defined by the three lowest mass possible 
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or probable cluster members by proper motion, BRB 2 7, 28 
and 29. Our estimate of a is in agreement with that from M03, 
0.60 + 0.11, and signifies that the number of objects per unit 
mass still increases at th ese low masses. In Fig. |7]the lognor- 
mal function obtained by Deacon & Hamblv ( 2004) is also rep- 
resented, scaled down so that its integral over the mass range 
0.5-0.026 Mq is equal to the number of objects in the survey 
area and in this mass range (56). 

Our estimate of a is similar to those for younger clusters 
as a Per (-80 Myr), IC 4665 (30-100 Myr) and A Orionis 
(~5 Myr). Their values are: 5 6 (0.2 > M(M Q ) > 0.06; 
Barrado v Navascues et all 12002). 0.6 (1 > M(M Q ) > 0.04; 
de Wit et alJl2006l) and 0.60 + 06 (0.86 > M(M Q ) > 0.024; 
Barrado v Navascues et alj E004). respectively. The presence 
of low-mass brown dwarfs after more than 100 Myr in the 
Pleiades open cluster, in numbers relative to stars which are 
similar to those found in other younger clusters, implies that 
differential evaporation of low-mass members relative to more 
massive ones has not been very significant - in contrast to 
what may occur in th e ~600 Myr old Praesepe open cluster 
JChappelle et alJEo05l) . Comparison of observations to evolu- 
tionary models of the Pleiades open cluster suggests indeed that 
only marginal different ial evaporation of t he massive brown 
dwarf s has occurred jMoraux et alJ 120041 iMoraux & C larke 
2005). Therefore searches for the lowest mass brown dwarfs 
and even planetary-mass objects should be conducted. A 
Pleiades giant planetary-mass object at the deuterium-burning 
mass limit has a predicted effective temperature of ~1300 K 
(DUSTY and COND models). This effectiv e temperature cor- 
respo nds to late L and early T field dwarfs JGolimowski et alJ 
2004) and magnitudes J > 20 ( Knap p et all2004l) when scaled 
to the distance of the cluster. Such a free-floating object in 
the Pleiades open cluster may be detected in a photometric 
search with mid-size telescopes (~4 m) but spectroscopic clas- 
sification will require a very large telescope of 8-10 m class. 
The discovery of such low-mass Pleiades objects will permit 
to compare with free-floating low-gravity candidates in the so- 
lar neighbourhood, as f or example the T6 brown dwarf SDSS 
Jl 11010.01+01 1613.1 tlKnappet alJl2004 . Extrapolating the 
mass spectrum slope to lower masses, we find that our sur- 
vey area may contain five brown dwarfs with masses 0.026- 
0.013 Mq and four planetary-mass objects with masses 0.013- 
0.05 Mq (or two per square degree). In the case where the mass 
spectrum is log normal down to planetary masses, these values 
drop to three and one, respectively (or one planetary-mass ob- 
ject per square degree). 

6. Conclusions 

As a result of an 77-band survey in the Pleiades, we have found 
18 brown dwarf candidates, most likely of L-type. They have 
17.4 < J < 19.7, / - J > 3.2 and theoretical masses -0.040- 
0.020 Mq. The near-IR follow-up for proper motion of 1 1 can- 
didates permits us to confirm six as cluster members. Another 
three remain as possible members. Two are clear non-members. 
The near-IR photometric follow-up of the seven other candi- 
dates indicates that at least five are probable contaminants, be- 
cause their J - K s colours are bluer than those of the cluster 



proper motion members. These latter determine the L brown 
dwarf photometric sequence of the Pleiades. The sequence in 
the J versus I - J colour-magnitude diagram at J > 18 is 
bluer than that predicted by the DUSTY atmospheres model, 
whereas it agrees with the predictions in the J versus J - H 
and J - K diagrams. Moreover the sequence overlaps that of 
the field L dwarfs shifted to the cluster distance in the three 
colour-magnitude diagrams, suggesting that the Pleiades and 
field L dwarfs may have comparable spectral energy distribu- 
tions and luminosities, and thus possibly similar radii. 

We find evidence for an intrinsic velocity dispersion of 
Pleiades brown dwarfs at least four times that of Pleiades stars 
more massive than the Sun. The estimated value >4 mas yr~' 
appears consistent with the kinematical expectations for brown 
dwarfs in a nearly relaxed cluster. 

Correcting for cont aminants and comparing w ith the proper 
motion objects from iDeacon & H amblv (2004) in the same 
fields, we obtain a substellar mass spectrum for the surveyed 
area of the cluster. For a power law dN/dM cc M~ a fit 
in the mass range 0.5-0.026 Mq, we find a spectral index 
a = 0.5 + 0.2 and in agreement with previous results. The slope 
is similar to that of much younger open clusters, suggesting the 
absence of significant differential evaporation of the low-mass 
brown dwarfs in the Pleiades. This supports deeper searches to 
detect lower mass free-floating objects in the cluster. 
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